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Abstract. Sharp bounds are obtained, under a variety of assumptions on the eigenvahies of 
the Einstein tensor, for the ratio of the Hawking mass to the areal radius in static, sphericaUy 
symmetric space-times. 

o 
o 

CN ; 1. Introduction 

> ! 

O . All of the space-times considered in this paper are connected, four-dimensional and satisfy 

w I the following conditions. 

^ , • Spherical Symmetry^ There is a time-like curve, called the time axis, with the property 

that at any point all normal directions are equivalent, i.e., that for any two space-like 
normal unit vectors there is an isometry of the space-time which fixes the point and 
takes the first vector to the second. This defines an action of 5*0(3) on the space-time 
bX)! whose orbits are called spheres. 

• Staticity. There is a one-parameter group of isometries, called time translations, whose 
CN ■ generating vector field is everywhere time-like. 

^ ■ • Regularity. The space-time, together with its metric, is of class C^, except possibly on 

'. 3-surfaces of discontinuity, where the second derivatives of the metric are allowed to 

^ I have jump discontinuities. 

The somewhat odd looking regularity assumption is borrowed from Lichnerowicz [9]. It is 
O ' meant to allow such discontinuities as one expects to find at the interface between two different 
^ ■ materials, but nothing worse. 

The areal radius r is defined by the requirement that the area of a sphere be Atct"^. In terms 
of the radius r and metric tensor g, we may then define the Hawking mass m by the relation 



cr: 



2m 



g^ikg^rdkr = 1 - ^ . (1.1) 

The purpose of this paper is to prove sharp bounds on the ratio 2m/ r under various hypotheses 
on the eigenvalues of the Einstein tensor. The particular hypotheses considered, their history 
and the resulting bounds are discussed in Section [31 

Three general comments should be made at this stage. First, the method employed is quite 
general and can be used to obtain sharp bounds on 2m /r for any matter model, not just those 
described below. Second, obtaining sharp bounds is, in each case, relatively easy. Proving 
sharpness, while not conceptually difficult, requires considerably more effort. Third, we carefully 
avoid the assumption, made tacitly by previous authors, that 2m < r. This point is discussed 
in more detail in the next section. 



""^We foUow Synge [lOj in caUing this assumption "spherical symmetry" for brevity. This is a bit misleading, 
as we are assuming more than just spherical symmetry. Our assumption excludes, for example, Schwarzschild 
space, which lacks a time axis. 
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Section [2] is devoted to a discussion of coordinates and the components of the Einstein tensor 
in our chosen coordinate system. Section [3] introduces the various assumptions on this tensor 
which are needed for the statement of our theorem. Our main resuh, Theorem 14.11 appears in 
Section HI while its proof is given in Section [5l 

2. Geometry and Coordinates 

A space-time of the class considered above has coordinates r, 6*, ^ and t, known as curvature 
coordinates, in which the metric takes the form 

gjk dxHx^ = e" dr^ + r\de^ + sin^ 6 d<^'^) - dt\ (2.1) 

Here, a and 7 are functions of r. As shown in Synge [10], the Einstein tensor in curvature 
coordinates is of the form 

= r-2_r-V"(l+r7'), 

G* = r-2-r-2e-"(l - to')- 

Here, primes denote derivatives with respect to r, while the off-diagonal entries are all zero. 
The formulae become a bit cleaner when one uses derivatives with respect to 

/3 = 2 log r, 

instead. Denoting such derivatives by dots, one obtains the equivalent system 



=r-2-r-V"(l + 27), (2.5) 
= = r-2e-"(-27 - f + d + «7), (2.6) 
G* = r-2 -r-V"(l -2d). (2.7) 

The corresponding Einstein tensor, given by Einstein's field equations, has diagonal entries 

Gl = -SirpR, Gl = Gl = -Sttpt, G\ = SiTfi (2.8) 



and all off-diagonal entries equal to zero. Here, p/j and pt are interpreted as the radial and 
tangential pressures, respectively, while fi is interpreted as the energy density. 
There are two annoying points about curvature coordinates. 

• The functions a and 7 may be of lower regularity than the metric, since r itself is of lower 
regularity than the metric. This is discussed in more detail by Israel [8j. For our purposes 
it suffices to note that regularity is, in the presence of the other assumptions, equivalent 
to the statement that a and 7 are functions of r, except possibly at certain points, 
where the radial pressure pr is continuous and the tangential pressure pt and energy 
density /i are allowed to have jump discontinuities. At r = 0, the correct condition is 
that a'(0) = y(0) = 0. 

• The coordinates may fail to cover the whole space-time. In fact, they cover the region 
from the time axis out to the first marginally trapped sphere, i.e., the first sphere where 
r = 2m. If we were to assume, as most authors do, that curvature coordinates cover 
the whole space-time, then we would, in effect, be making the very strong additional 



(2.2) 
(2.3) 
(2.4) 
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assumption that 2m /r < 1 everywhere. This we wish to avoid. For the classes of space- 
times we consider it is, in fact, the case that 2m /r < 1 everywhere, but this belongs to 
the conclusion of our theorem, not to its hypotheses. 

The simplest example of a space-time that satisfies our spherical symmetry, staticity 
and regularity assumptions but has a marginally trapped surface is de Sitter space, for 
which e" = —e~^ = 1 — jW' . In this case, the coordinates cover a region where r < R 
but break down at the boundary. Outside this region, there is another which is isometric 
to the first, and it is easy to check that 2m = r at r = R. However, de Sitter space does 
not satisfy the hypotheses of our theorem because it has negative pressures everywhere. 

3. Matter Models 

Various conditions on the three functions Pr, Pt and /i are of interest: 

• Non-negative Isotropic Pressure: For fluids, one expects Pr = pr ^ 0. The sharp bound 

2m/r < I2V2 - 16 ^ 0.9706 

under this assumption, and no others, was derived by Bondi [1]. His method of proof is 
closely related to ours but is not rigorous. 

• Buchdahl Assumption: For static stars with constant density, one has the bound 

2m/r < 8/9 

derived by Buchdahl [5j. More generally, this bound holds when p_r = Pt ^ 0, as long 
as /i > is decreasing; see |i5j. The isotropy assumption was relaxed in [7J, where the 
case J9r > > was treated, still the monotonicity assumption remains crucial. 

• Dominant Energy Condition: For almost any reasonable matter model, one expects 
\Vr\ ^ S'lid \pt\ < /U- In the special case that Pr,Pt > 0, the bound 

2m/r < 48/49 ^ 0.9796 

is provided by Our bound for this special case is roughly 0.963, which we show to 
be sharp. 

• Vlasov- Einstein: For Vlasov-Einstein matter, the stress energy tensor is an integral of 
those of individual particles, each of which has rank one and satisfies the dominant 
energy condition. This implies that Pr > 0, pt > and pr + 2pT < yU. Under these 
assumptions, Andreasson [5] has recently shown that the sharp@ bound is 

2m/ r < 8/9. 

Our method provides a new, and considerably shorter, proof of this result. 

• Zero Radial Pressure: The case Pt > P_r = was studied by Florides ^ who obtained 
the sharp bound 

2m/ r < 2/3. 

This can also be proved using our method, but the resulting proof is neither shorter nor 
clearer than the original, so we do not consider this case further. 

somewhat unfortunate feature of our argument is that the sharpness of the estimate ^ < | is proved 
only within the class of space-times satisfying the pressure conditions above, without considering whether such 
space-times arise from solutions of the full Vlasov-Einstein system. Andreasson's argument, on the other hand, 
does provide solutions to the full system. 
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4. Our main result 

Theorem 4.1. Consider a space-time satisfying the regularity, staticity and spherical symmetry 
conditions described in the introduction. Suppose that the corresponding Hawking mass fll.ll) is 
finite and that the pressures Pr,Pt and energy density /i are all non-negative. 

(1) Vlasov- Einstein case. Assuming that p^ + 2pT < fi, one has 

6m o \^ f 2m\ 2m 



A- — + 87rr^Pnj < 16 ^1 - — J , — <-. (4.1) 

(2) Isotropic case. Assuming that p^ = px, one has 

( 2m 9 \^ / 2m \ 2m r- 

( 2 - — + Svrr^p^ j < 36 f 1 - — j , —< 12^2 -16. (4.2) 

(3) Isotropic case with dominant energy. Assuming that p a = pr ^ fJ', one has 

(2- — + 87rrO ' < 9.551 (l-—), —< 0.865. (4.3) 
\ r J V ^ / 

(4) Dominant energy in tangential direction. Assuming that px < /i, one has 

10m . \^ / 2m\ 2m 2^2 + 2 , , 

6 + ST^r^pR < 40 1 , — < . (4.4) 



r / V / ^ 5 

(5) Dominant energy case. Assuming that p^jp^ < fi, one has 

10m 9 \^ / 2m\ 2m 

6 + STcr^pR < 37.924 1 , < 0.963. 4.5 

r / \ r J r 

Moreover, these ten estimates are all sharp for the class of space-times considered in each case. 
As for the numerical values that appear in (14.31) and (14. 5p . these can be described in terms of a 
system of ODEs which arises in the course of the proof; see (15.241) . The values given here are 
accurate up to three decimal places. 

Remark 4.2. There is no assumption on the behavior of the space-time as r tends to infinity. 
In fact, we do not even assume that r is unbounded. This point is crucial. It allows us to 
apply the theorem to the interior of a finite sphere and, in particular, to the interior of the first 
marginally trapped surface, if there is such a surface. 

More precisely, suppose we can prove the theorem in the region where the curvature coordi- 
nates are defined, namely in the region where 2m < r. For each matter model, we may then 
deduce that 2m /r < c < 1 for some constant c which depends on the matter model considered. 
Since 2m /r is continuous and our space-time is connected, this actually implies that 2m /r < c 
throughout the space-time. In other words, the marginally trapped surface that we allowed 
is not, in fact present, and the curvature coordinates, which might a priori have covered only 
part of the space-time, cover the whole space-time. We therefore obtain the full theorem from 
the special case where the whole space-time is covered by curvature coordinates. In particular, 
we may, and do, use curvature coordinates throughout the proof without further comment. 

Remark 4.3. Our proof for the isotropic case pr = pt applies verbatim in the more general 
case pr > Pt, while the estimates (14.21) are sharp for that case as well. 
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Remark 4.4. The assumptions of Theorem 14.11 can be shghtly improved in the sense that we 
do not use our hypothesis px > to estabhsh the given estimates. This hypothesis is merely 
included to improve the conclusions of Theorem 14.11 as sharpness is now shown over a smaller 
class of space-times. In fact, the examples we construct in order to prove sharpness belong to 
the even smaller class of space-times which are vacuum outside a sphere. 

The proof of Theorem 14. II is based on the following elementary fact, which is essentially due 
to Bondi [1]. 



Lemma 4.5. Let the assumptions of Theorem 4_J_ hold. Then the variables 

Off) 

x = 1 - e-° = — , y = -r^Gl = Snr^pn (4.6) 

give rise to a parametric curve which lies in [0, 1) x [0, oo) and satisfies the equations 

Snr^PR = y, (4.7) 

'"^^-= 2(1^^ + ^+4(1^' ^'-'^ 
SvrrV = 2x + x, (4.9) 
where the dots denote derivatives with respect to (3 = 21ogr. 

Proof. First of all, we combine equations (12. 4p and (12. 8p to write 

87rrV = l-<9^(re~"). 
Integrating over [0, r] and using the definition of the Hawking mass (11.11) . we then get 

— = l-e-" = x. (4.10) 
r 

This implies x > because m > whenever /x > 0. Next, we use (12. Sp to get 

y = -r^G; = Srrr^pR > 0. 
To establish our assertion (14. 9p . we combine (12. Sp . (12. 7p and (I4.10p to find that 

SvrrV = r^Gl = 1 - e~"(l - 2d) = 1 - (1 - x) (^1 - =x + 2x. 

To establish our remaining assertion (14. 8p . we first use (12. 5p and (I4.10p to get 



y = -r'G: = -l + {l-x){l + 2'j), a-- 

1 — X 

Solving the leftmost equation for 7 and differentiating, we conclude that 

x + y {1 + y)x + {1 - x)y 

^'^ = 1 ' ^'^ = n \2 ■ 

1 — a; (1 — xj^ 

On the other hand, equations (12.80 . (12.60 and (I4.10p combine to give 

Snr'^PT = (1 — x)(27 + 7^ — d — 0:7). 

Using these facts and a simple computation, one may thus easily deduce (14. 8p . 
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5. Proof of Theorem 14.11 

To prove the desired estimates, we study the curve (14.61) provided by Lemma 14.51 In each 
case, we are seeking an upper bound for x = 2m /r and also an upper bound for 

(n(l - x) + l + yy 



1 



where the exact value of n varies from case to case. Differentiating (15. ip . we get 

n(l — x) + 1 + y 



:i-x) 



[1 + y — n{l — x))x + 2(1 — x)y 



(5.1) 



(5.2) 



throughout the curve (14. 6p . where dots denote derivatives with respect io (3 = 21ogr. In the 
special case that n = 1, this formula reads 

2 — X + y 



Wi 



{x + y)x + 2(1 — x)y 



(5.3) 



81). Let us also recall that 



(1-X)2 

and it is closely related to the tangential pressure pt', see 

0<a;<l, y>0 
throughout the curve (14.60 . a fact we shall frequently need to use in what follows. 

5.1. Vlasov-Einstein case. In this case, we are assuming that pn + 2pT < /U. According to 
Lemma [4. 5 [ the corresponding curve (14.60 must thus satisfy 

Z3ix,y) 



{3x + y - 2)x + 2{1 ~ x)y < 



Z3 



3x^ -2x + y^ + 2y. 



(5.4) 



Combining the last equation with our computation (15.21) . we now find 

4 — 3x + y 



W3 



X] 



{3x + y - 2)x + 2{l - x)y 



A-3x + y 
- 2(l-a;)2 '^^ 

In particular, is decreasing whenever ^3 > 0, so it must be the case that 

W3 < max wz{x,y) = ^3(0,0) = 16 

0<x<l 
Z3<0<y 

throughout the curve. This proves the first inequality in (14.11) . which also implies the second 
inequality because the maximum value of x over the region 0<x<l,?/>0, W3<16is 
attained at (8/9,0), namely at the point at which the curve W3 = 16 intersects the x-axis. We 
refer the reader to Fig. 1 for a sketch of the curves = and = 16. 

To show that the estimates in (14. ip are sharp, we need to construct a space-time such that 
the corresponding curve of Lemma intersects a small neighbourhood of (8/9, 0). Let us now 
temporarily assume that we have a parametric curve 

x = x{t), y = y{T); r G (0, oo) 

which passes near the point (8/9,0) and also satisfies the following properties: 
(Al) J- ■ is both negative and integrable; 
(A2) 0< x\t) < 1 and ?/(r) > for each r > 0; 
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(A3) y(r) = for all large enough r and x(r) — > as r — > oo; 
(A4) the curve is except for finitely many points. 

Given such a curve, we can easily construct a space-time as follows. First, we define 

^ dws/dr 2{l-xy 

^ ' zz{x,y) A-?,x + y ^ - ' 

and we note that k is both positive and integrable by (A1)-(A2). Next, we define 

13 = j KdT, r = exp(/?/2) (5.6) 

and finally, we define the metric coefficients in (12. ip by 

/X -\- y 

^^Y— -y-Krfr. (5.7) 

Letting dots denote derivatives with respect to [3 = 21ogr, as usual, we then get 

1 dw^ 4 — 3x + y 
K, dr 2 1 — xr 



using our definitions (15. 6p and (15.51) . In view of our computation (15. 2p . this gives 



{3x + y- 2)x + 2(1 - x)y = -'J^ , (5.8) 

which is equivalent to the equation + 2^^ = n because of Lemma 14.51 

To finish the proof for this case, it thus remains to construct the curve whose existence we 
assumed in the previous paragraph. We have to ensure that the curve satisfies (Al)-(A4), that 
it passes arbitrarily close to (8/9,0) and that the corresponding quantities p_r, Pt,/^ provided 
by Lemma 14.51 are non-negative. Let us then fix some small e > and consider the curve 

W3-3e(x,?/)= eVl + 3a; + 4(l -£) ^ (5.9) 

When £ = 0, this reduces to the curve ^3 = 16 which passes through the origin and (8/9,0). 
When e = 1, it reduces to the curve Z3 = which passes through the origin and (2/3,0). In 
the more general case < £ < 1, it describes a curve that lies between these two curves. We 
start out at the origin and follow this curve until we hit the x-axis, and then we return to the 
origin along the x-axis. Let us henceforth denote by Ci the curve obtained in this manner; we 
refer the reader to Fig. 1 for a typical sketch of this curve. 

The fact that Ci satisfies (A2)-(A4) is trivial. To check that it satisfies (Al) along the part 
defined by (15. 9p . we recall that this part lies between the curves 1^3 = 16 and ^3 = 0. Thus, it 
is easy to see that 2:3 > along this part, and we need only check that 

^ < (5.10) 

dr 

as one follows the curve (15.91) in the positive x-direction. Differentiation of (15. 9p gives 

dW3_3s _ 3£:^W3-3g dx 

dr y/1 + 3x dr 

3e 3(1 -£)(1 -x) + 1 + ?/ dx 
VI + 3x VI - X dr 
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Figure 1. The curve Ci for the Vlasov-Einstein case. 



along the curve fl5.9p . and we may compare this equation with (15. 2 p to find that 



2(1-0;) 



dr 



3g(l-a;)3/2 



3(1-£)(1 



X) 



1-y 



dx 
dr 



(5.11) 



along the curve (15. 9p . Employing our computation (15. 2p once again, we deduce that 

dw^ 



3£(4 - 3x + y) 



VT+~3x dx 
dr 



dr I- X VI + 3a; 

Since ^ > here, the desired (15.100 follows. To show that (Al) also holds for the remaining 
part of the curve Ci, we need only note that 

1 dw?j 1 4 — 3x 



dx 



dr X {1 — X) 

along the line y = because this line is traversed in the direction of decreasing x. 

Finally, we check that PR,PT,fJ' > throughout the curve Ci. The fact that Pr> follows 
by (A2) because STir^pR = y hj definition. Since (15.80 ensures that fi = Pr + 2pT, we need only 
check that Pt > as well. Let us now write 

{x + yf 



It^t'^Pt 



X + y 
2(1 -x) 
1 — a; 



X + y + 



Wi + 



4(1 -x) 

{x + yf 



2(2- x + y) 4(1 -a;) 
using equations (14. 8 p and (15. 3p . Along the part of Ci defined by (15.90 . we have 



(5.12) 



Wi 



2 — X + y 
1 — X 



X 



+ 2-3e 



X 



by (15. Sp and (15. lip . In view of our definition (15. 6p . we thus have 



Wi 



2 — X + y 
1 — X 



X 



. Vl + 3x 



+ 2 - 3e 



dx/dr 



K 
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Since e > is small and k is positive by above, this implies Wi > 0, hence px > hy (I5.12p . 
For the remaining part of Ci along the x-axis, Lemma 14.51 and (15. 8p give 

XLi 

PR = 0, pt = r , /i = 2pT, 

4(1 — X) 

so it easily follows that Pr = pr = fJ' = throughout this part of the curve. 
5.2. Isotropic case. In this case, our assumption that = px is equivalent to 

(x + i/)x + 2(l-x)2/ = -^i^y^, zi = {x + yf -Ay{l-x). (5.13) 
Proceeding as before, we use our computation (15. 3p to find that 



2 — X + y 



{x + y)x + 2{l - x)y =-777- -^■zi{x,y). (5.14) 

2 1 — xY 



Once again, wi is decreasing as soon as 2:1 > 0, so it must be the case that 

Wi < max wi{x,y) = Wi(0,4) = 36 

0<x<l 
zi<0<y 

throughout the curve. This proves the first inequality in (14.20 . while the second inequality 
follows because the maximum value of x over the region 0<x<l,?/>0, Wi<36is attained 
at (12x/2- 16,0). 

To show that the estimates in (14. 2 p are sharp, we argue as in the previous case. Suppose we 
have a curve which passes near the point (12v^ — 16, 0) and satisfies 

(Bl) ^ ■ is both negative and integrable 

as well as (A2)-(A4). Then we can follow our previous approach with 

.(r) = . ^ii^ > (5.15) 

zi{x,y) 2-x + y 

instead of (15. 5p . Our definitions (I5.6p - (l5.7p are still applicable, however they now imply 

1 dwi 2 — X + y 
K, dr 2(1 — X 



In view of our computation (15. 3p . they thus imply 

(x + y)x + 2(1 - x)y = -£ii|l^ , (5.17) 

which is equivalent to the equation p^ = p^ because of Lemma 14.51 

To finish the proof for this case, it thus remains to construct the curve whose existence we 
assumed in the previous paragraph. Fix some small e > and set 

Xs = e, ye = 2~3e + 2y/{l-e){l-2e) (5.18) 

for convenience. Then [x^^ye) is a point on the curve = which is close to (0,4). To define 
the first part of the desired curve, we use the equation 

Vw^x^ = V^m^IM^ + 2^a;,x - A^x^ (5.19) 
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y 




Figure 2. The curve C2 for the isotropic case. 



where is determined by requiring that the curve passes through {xe,ye), namely 



A^= ^Mys)-2 _ ^^^^^^ 



e 



We start out at the origin and we fohow the curve fl5.19p until we reach the point (x^, y^), then 
we follow the curve 

/ / six — X 

\/wi{x,y) = ^Jwi{xe,ye) , x>Xs (5.21) 

V 1 - a; 

until we hit the x-axis, and finally we return to the origin along the x-axis. Let C2 denote the 
curve obtained in this manner; a typical sketch of this curve appears in Fig. 2. 

The fact that C2 satisfies (A2)-(A4) is trivial; we now check that it satisfies (Bl). When it 
comes to the part of C2 defined by f l5.19p . we have zi < 0, x < Xe and 

dwi/dr , , ^ dx , ^ 

— ^ = 2A, x,-x ■ — . 5.22 
2y/wi dr 



Since x is increasing along this part of C2, it thus suffices to check that is positive. In view 
of fl5.20p . this is certainly the case for all small enough e > because 



lim e^Ae = ^^1(0, 4) -2 = 4. 



When it comes to the part of C2 defined by f l5.2ip . we have Zi > 0, x > x^ and 

dwi/dr , , 4(1 — x) + X — Xe dx 

= —e x — xA ■ — — -, — - — TTTTT, ■ -;- < 0, 

2^^ ^ ^' 2(1 -x)^/2 dT - ' 

as needed. When it comes to the remaining part of C2 along the x-axis, we have 

dw\ /dr 2 — X dx ^ 
Zi x(l — xY dr ' 

and this shows that the desired property (Bl) holds throughout the curve €2- 
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Finally, we check that Pr,Pt,I^ > throughout the curve C2. The fact that Pr > follows 
trivially as before, hence = P_r > by fl5.17p and we need only worry about /i. Since 

Svrr^/i = 2x + X 

by (14.91) . we have /i > as long as x is increasing along the curve, so we need only check the 
part of C2 along the x-axis. As in the previous case, however. Lemma 14.51 and fl5.17l) combine 
to give Pr = Pt = fJ' = throughout this part, so the proof for this case is complete. 

5.3. Isotropic case with dominant energy. Our assumption that Pij = Pt < A* gives 

{x + y)x + 2{1 - x)y = -^^^^Y^ , 2x > y - x (5.23) 

with Zi as in (I5.13p . Due to the isotropy condition, (15.141) remains valid, so Wi is increasing 
if and only if Zi < 0. Since the curve of Lemma [4.51 starts out at the origin, where Zi = 0, it 
may only attain the largest possible value of wi at a point along the curve zi = 0. It is easy 
to check that higher values of wi occur at higher points on this curve. To attain the largest 
possible value of Wi, the curve of Lemma 1431 must thus ascend as fast as possible within the 
region Zi < 0. Since it starts out at the origin, it must satisfy 

Zi (x, I/) 

(x + y)x + 2(1 — x)y = ^ — , 2x = y — x 

until it exits the region zi < 0. This gives rise to the system of ODEs 

2x = y-x, 2y= ~ ~ (5.24) 

1 — X 

which has a saddle point at the origin. The solution of interest is the one corresponding to 
the unstable manifold associated with the origin. Using numerical integration, we find that it 
intersects the curve 2:1 = at the point {xi,yi) = (0.4927,0.6939); see Fig. 3. This makes 

Wi{xi,yi) ~ 9.551 

the largest possible value of Wi, and then we can use the fact that Wi < 9.551 to deduce that 
the largest possible value of x is attained at (0.865, 0). 

To show that our results for this case are sharp, we need to find a curve which passes near 
the point (0.865,0) and satisfies (Bl) as well as (A2)-(A4). Given such a curve, one can use 
our approach in the previous case to obtain a space-time for which Pr = px- We start out at 
the origin and we follow the solution to the ODE 

dy ^ y{2 - 3x - y) ^5) 
dx (1 — x){y — x) 

corresponding to the associated unstable manifold; we do so until we reach the point (xi,?/i) 
that lies on the curve zi = 0, then we follow the curve 

^ywl{x,y) = a/wi(xi,|/i) , x>xi (5.26) 

V 1 — a; 

until we hit the x-axis, and finally we return to the origin along the x-axis. We refer the reader 
to Fig. 3 for a typical sketch of the curve C3 obtained in this manner. 
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Figure 3. The curve C3 for the isotropic case with dominant energy. 

The only nontrivial properties we need to verify are (Bl) and the fact that < /x. When it 
comes to the part of C3 defined by f l5.25p . we have Pr = Pt = and also 

i_ . ^ = 2-x + y . ^ < g 27) 

Zi dr (1 — — x) dr ~ ' 

so the desired properties are easily seen to hold. The same is true for the part of C3 along the 
X-axis because Pi? = Pt = = and since 

1 dwi 1 2 — X <^ Q 

Zi dr X (1 — x)^ dr 



along this part. For the remaining part defined by fl5.26p . we have 

dwi/dr 4(1 — x) + X — xi dx 



—E[X — XiJ ■ 



< 0, 



20nr 2(l-x)3/2 dr 

which implies property (Bl) because Zi > for this part. Writing 05.161) in the form 



(5.2^ 



Wi 



2-x + y ^/wi{x,y) 

zi{x, y) = r^TT^ ■ zi[x, y), 



2(1 -x)2 -^^-^'^^ 2(l-x)3/2 
we now combine the last two equations to deduce that 

e{x — xi)(4 — 3x — xi) 

throughout the curve 05.26P . According to Lemma 14. 5[ the condition Pr < n we need to verify 
is equivalent to the condition 2x > y — x, so we need to check that 

zi{x,y) 



> e{A — 3x — Xi){y — x) 

X — Xi 

throughout the curve 05.26p . Write equation 05.260 in the equivalent form 
y = /(x) = A/ti?i(xi,?/i)Vl - X + X - 2 - e(x - xi)^ 



(5.29) 
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Then Zi{xi, /(xi)) = Zi{xi,yi) = by construction, so one easily finds 

Z1{XJ{X)) o/ , N . /g , 0\ _ . 7^fi 

lim = 8(xi + — 4 — (3xi + Hi — 2) ■ — — ^ 4.746 

x^xi X — Xi ^/l — Xi 

using fl5.13p . Thus, the left hand side of fl5.29p is bounded away from zero near x = Xi. Since 
the same is true away from x = Xi, where zi itself is bounded away from zero, we can always 
find a small enough e > so that (15.291) holds throughout the curve fl5.26p . 



5.4. Dominant energy in tangential direction. Our assumption that Pt < gives 

(5x + ?/-4)x + 2(l-x)?/ < -^^^y^, z^ = {x + yf -4:x{l-x). (5.30) 
Proceeding as before, we use our computation (15. 2p to find that 



Q — hx + y 

W5 = 



i5x + y- 4)x + 2(1 - x)y < -^—^^±1 ■ z,ix, y). (5.31) 

2(1 — xY 



1 — X 

Once again, is decreasing as soon as 2:5 > 0, so it must be the case that 

< max W5(x, y) = W5(l/10, 1/2) = 40 

0<a;<l 
Z5<0<y 

throughout the curve. This proves the first inequality in (14.40 . and the second inequality follows 
as before. 

To show that the estimates in (14.40 are sharp, we need to find a curve which satisfies 

(CI) ^ ■ ^jp- is both negative and integrable 

as well as (A2)-(A4). Given such a curve, one can use our previous approach to obtain a 
space-time for which px = fJ^. To define the first part of the curve, we use the equation 

Ax 



^/w5ix, y) = v/«^5(0, 0) + — - Ax\ (5.32) 
where A is chosen so that the curve passes through (1/10, 1/2), namely 

A= 100(\/40-6) > 0. 

We start out at the origin and we follow the curve (15.321) until we reach the point (1/10, 1/2), 
then we follow the curve 

^Jw^ix, y) = ^7/75(1/10, 1/2) - ^^"^ -^Q)' ^ x>l/10 

Vl - X 

until we hit the x-axis, and finally we return to the origin along the x-axis. Since this curve is 
almost identical with the one for the isotropic case, our previous approach applies with minor 
changes; we shall not bother to include the details here. 



14 



PASCHALIS KARAGEORGIS AND JOHN G. STALKER 



5.5. Dominant energy case. In this case, our assumption that Pr,Pt ^ 1^ gives 

{5x + y - A)x + 2{1 - x)y < -^^^^Y^ , 2x > y - x (5.33) 

with as in f l5.30p . Since f l5.3ip remains vahd, is decreasing when > 0, so its maximum 
value is attained in the region 2:5 < 0. To obtain the largest possible value of w^, we need to 
ensure that is as large as possible in this region. In view of (15.311) . this simply means that 
equality must hold in the first inequality in fl5.33p . We are thus faced with a situation which is 
almost identical with fl5.23p . Arguing as before, we find that the curve must satisfy 

2x = y — X, (5x + y — 4:)x + 2(1 — x)y = — ^^2^^ 

until it exits the region z^ < 0. This is the same system of ODEs that we had in fl5.24p . and 
the rest of our argument applies almost verbatim. The solution associated with the unstable 
manifold at the origin intersects the curve = at the point (0.2746, 0.6180) and so 

^5(0.2746, 0.6180) ^ 37.924 

is the largest possible value of 1^5. Using this fact, we get the upper bound on x which is stated 
in the theorem. To show that our results for this case are sharp, we follow our approach in the 
isotropic case with dominant energy. As there are only minor changes that need to be made, 
we are going to omit the details. 
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